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ABSTRACT 

A model is derived for relativistic Foynting jets from the inner region of a disk around a rotating 
black hole which is initially threaded by a dipole-like magnetic field. The model is derived from the 
special relativistic equation for a force-free electromagnetic field. The "head" of the Foynting jet is 
found to propagate outward with a velocity which may be relativistic. The Lorentz factor of the head is 
F ~ [Bq/ {SiTrTZ'^ pextc^]^^^ if this quantity is much larger than unity. For conditions pertinent to an active 
galactic nuclei, F w 8{10/ny/^{Bo/10^Gy/^{l/cm^ /uexty^^ , where Bq is the magnetic field strength 
close to the black hole, pext = fhriext is the mass density of the ambient medium into which the jet 
propagates, TZ = r^/rg > 1, where is the gravitational radius of the black hole, and rg is the radius of 
the O— point of the initial dipole field. This model offers an explanation for the observed Lorentz factors 
~ 10 of parsec-scale radio jets measured with very long baseline interferometry. 

Subject headings: galaxies: nuclei — galaxies: magnetic fields — galaxies: jets — quasars: general 
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1. INTRODUCTION 

Highly-coUimated, oppositely directed jets are observed 
in active galaxies and quasars (see for example Bridle 
& Eilek 1984; Zensus, Taylor, & Wrobel 1998), and in 
old compact stars in binaries (Mirabel & Rodriguez 1994; 
Eikenberry et al. 1998). Further, well coUimated emission 
line jets are seen in young stellar objects (Mundt 1985; 
Biihrke, Mundt, & Ray 1988). Recent work favors models 
where the twisting of an ordered magnetic field threading 
an accretion disk acts to magnetically accelerate the jets 
(e.g., Meier, Koide, & Uchida 2001; Bisnovatyi-Kogan & 
Lovelace 2001). There are two regimes: (1) the hydromag- 
netic regime, where energy and angular momentum are 
carried by both the electromagnetic field and the kinetic 
flux of matter, which is relevant to the jets from young 
stellar objects; and (2) the Poynting flux regime, where en- 
ergy and angular from the disk are carried predominantly 
by the electromagnetic field, which is relevant to extra- 
galactic and microquasar jets, and possibly to gamma ray 
burst sources. 

Different theoretical models have been proposed for 
magnetically dominated or Foynting outflows (Newman, 
Newman, & Lovelace 1992) and jets (Lynden-Bell 1996, 
2003) from accretion disks threaded by a dipole-like mag- 
netic fleld. Later, stationary Foynting flux dominated 
outflows were found in axisymmetric magnetohydrody- 
namic (MHD) simulations of the opening of magnetic loops 
threading a Keplerian disk (Romanova et al. 1998). MHD 
simulations by Ustyugova et al. (2000) found coUimated 
Foynting flux jets. The present work represents a con- 
tinuation of the studies by Li et al. (2001) and Lovelace 
et al. (2002; hereafter L02), which are closely related to 
the work by Lynden-Bell (1996). Self-consistent force-free 
fleld solutions are obtained where the twist of each field 
line is that due to the differential rotation of a Keplerian 
disk. 



2. THEORY OF POYNTING OUTFLOWS 

Here, we consider relativistic Foynting outflows from a 
rotating accretion disk around a Kerr black hole. We as- 
sume that at an initial time f = an axisymmetric, dipole- 
like magnetic field threads the disk. The initial field ge- 
ometry is shown in the lower part of Figure 1. This field 
could result from dynamo processes in the disk (e.g., Col- 
gate, Li, & Pariev 1998). A more general magnetic field 
threading the disk would consist of multiple loops going 
from radii say Va to rb > ra and Tc > rf, to > Tc, 
etc. It is not clear how to treat this case analytically, but 
non-relativistic MHD simulations of this type of field show 
strong Foynting fiux outfiows (Romanova et al. 1998). 

The disk is assumed to be highly conducting and dense 
in the sense that pd{r^)^ ^ B^/47r, pd is the density and 
ri(r) the angular rotation rate of the disk. Thus a magnetic 
field threading the disk is frozen into the disk. Further, we 
suppose that the radial accretion speed of the disk is much 
smaller than azimuthal velocity fir. For a corotating disk 
around a Kerr black hole 



(1) 



for r > Tms where rms is the innermost stable circular or- 
bit, where a* is the spin parameter of the black hole with 
> a* < 1 and rg = GM/c?. For example, for a* = 0.99, 
fms ~ 1.45rg. The simplification is made that equation 



(1) also appfies for r < 



The region r < 



has a 



negligible infiuence on our results for the considered con- 
ditions where the radial scale of the magnetic field tq is 
such that {rg/r^Y <^ 1. 

In the space above (and below) the disk we assume 
a "coronal" or "force-free" ideal plasma (Gold & Hoyle 
1960). This plasma may be relativistic with fiow speed v 
comparable to the speed of light. Away from the head 
of the jet, the electromagnetic field is quasi-stationary. 
This limit is applicable under conditions where the energy- 
density of the plasma 7pc^ is much less than the electro- 
magnetic field energy-density (E^ -I- 'B'^)/8n. 
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Fig. 1. — Sketch of the magnetic field configuration of a Poynting jet. The bottom part of the figure shows the initial dipole-like magnetic 
field threading the disk which rotates at the angular rate fl{r). The top part of the figure shows the jet at some time later when the head of 
the jet is at a distance Z{t). At the head of the jet there is force balance between electromagnetic stress of the jet and the ram pressure of the 
ambient medium of density pext ■ 



Cylindrical (r, </>, z) coordinates are used and axisym- 
metry is assumed Thus the magnetic field has the form 
B — Bp + B^(p , with Bp — Bj-y + B^i ■ We can 
write Br = -{l/r){d^ /dz), B^ = (l/r)(5*/ar), where 
\['(r, z) = rA^{r,z) is the flux function. In the plane of 
the disk, the flux function is independent of time owing the 
frozen-in condition. A representative form of this function 
is 
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vl/(r,0) = - 



r'Bo 



2 1 + 2(r/ro)3 
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where Bq is the axial magnetic field strength in the cen- 
ter of the disk, and tq is the radius of the O— point of 
the magnetic field in the plane of the disk as indicated in 
Figure 1. Equation (2) is taken to apply for r > even 
though it is not valid near the horizon of the black hole. 
As already mentioned the contribution from region is neg- 
ligible for the considered conditions where [rg/rof' <C 1. 
Note that ^(r, 0) has a maximum, with value r^Bo/G, 
at the O-point where B(r'o,0) = 0. For (r/ro)^ < 1, 
^'(r, 0) (X which corresponds to a uniform vertical field, 
whereas for (r/ro)^ ^ 1, ^(r, 0) oc l/r which corresponds 
to the dipole field B^ir, 0) oc -l/r^. 

It is clear that equations (1) and (2) can be combined 
to give n — ri(^') which is a double valued function of 
5*. The upper branch of the function is for the inner part 
of the disk (r < tq), while the lower part is for the outer 
part of the disk. A good approximation for the upper 
branch is obtained by taking « r^Bo/2 which gives 



n = (cVGAf)/[a, + (2vI//i?or2)3/4]. 

The main equations for the plasma follow from the con- 
tinuity equation V-(pv) = 0, Ampere's law, V x B = 
47rJ/c, Coulomb's law V • E = 47rpe, with pe the charge 
density, Faraday's law, V x E — 0, perfect conductiv- 
ity, E-|-vxB/c = 0, with v the plasma flow veloc- 
ity, and the "force-free" condition in the Euler equation, 
PeE -|-JxB/c = (see for example Lovelace, Wang, 
& Sulkanen 1987). Owing to the assumed axisymmetry, 
= 0, so that the poloidal velocity Vp = kB^. Mass 
conservation then gives B • V(pk) — 0, which implies that 
pK — F(^')/47r, where F is an arbitrary function of 5". In 
a similar way one finds that — kB^ = rG'(^), so that 
E = -G(*)V4', and rB^ = so that there are two 

additional functions, G and H. 

The function G is determined along all of the field lines 
which go through the disk. This follows from the perfect 
conductivity condition at the surface of the disk, z = 0, 
Er+{-v X B)r/c = 0. This gives = —{v^Bz—VzB^)/c = 
—v^Bz/c, where Vz is zero at the disk and is the disk 
velocity. Therefore, Er{r,0) = -D, [d*(r, 0)/dr]/c, so that 
G{^) = n{r)/c which gives n = 

The component of the Euler equation in the direction of 
gives the force-free Grad-Shafranov equation, 



with A* 



2r^ 



HH' = , (3) 



d^/dr^ - {l/r){d/dr) + d^/dz^ and H' 
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dH {'i') / d'i' (Scharlcmann & Wagoner 1973; Lovelace, 
Wang, & Sulkanen 1987). 

As mentioned above, we consider an initial value prob- 
lem where the disk at t = is threaded by a dipole-Hke 
poloidal magnetic field. The form of H{'i>) in equation (3) 
is then determined by the differential rotation of the disk: 
The azimiithal f,wist of a given field line going from an in- 
ner footpoint at ri to an outer footpoint at r2 is fixed by 
the differential rotation of the disk. 

Ampere's law gives ^rfl • B = (Att/c) J dS ■ J, so that 
rB^{r,z) = H{'^) is (2/c) times the current flowing 
through a circular area of radius r (with normal z) labeled 
by ^'(r, z)= const. Equivalently, —H[^{r, 0)] is (2/c) times 
the current flowing into the area of the disk < r. Our pre- 
vious work (Li et al. 2000; L02) shows that -H{^') has 
a maximum so that the total current flowing into the disk 
for r < r„ is lioL = {c/2){-H)max, where is such that 
—H[^{rm, 0)] = {—H)max where is less than the radius 
of the O— point, vq. The same total current Itot flows out 
of the region of the disk r = r„i to r^. 

For a given field line we have rd(j)/B^ = dsp/Bp, where 
dsp = V dr^ + dz^ is the poloidal arc length along the field 
line, and Bp = ^^B^ + B^. The total twist of a field line 
loop is 




with the sign included to give A(/) > 0. The integra- 
tion goes from the disk at a radius ri < tq out into 
the corona and back to the disk at a radius r2 > tq. 
For a prograde disk the field line twist after a time t is 
A<t>{-^) = n{ro)t{[n{ri)/n{ro)] - [r!(r2)/f^(ro)]}. 

2.1. Poynting Jets 

Our previous study of non-relativstic Poynting jets by 
analytic theory (L02) and axisymmetric MHD simulations 
(Ustyugova et al. 2000) showed that as the twist, as mea- 
sured by fl{ro)t, increases a new, high twist field config- 
uration appears with a different topology. A "plasmoid" 
consisting of toroidal flux detaches from the disk and prop- 
agates outward. The plasmoid is bounded by a poloidal 
field line which has an X— point above the O— point on 
the disk. The occurrence of the X— point requires that 
there be at least a small amount of dissipation in the evo- 
lution from the poloidal dipole field and the Poynting jet 
configuration. The high-twist configuration consists of a 
region near the axis which is magnetically collimated by 
the toroidal field and a region far from the axis which 
is anti- collimated in the sense that it is pushed away from 
the axis. The field lines returning to the disk at r > rg 
are anti-coUimated by the pressure of the toroidal mag- 
netic field. The poloidal field fills only a small part of the 
coronal space. 

In the case of relativistic Poynting jets we hypothesize 
that the magnetic field configuration is similar that in the 
non-relativistic limit (L02, Ustyugova et al. 2000). Thus, 
most of the twist Acj) of a field line of the relativistic Poynt- 
ing jet occurs along the jet from z = to Z(t) as sketched 
in Figure 1, where Z{t) is the axial location of the "head" 
of the jet. Along most of the distance z = — Z the 
radius of the jet is a constant and ^ = ^{r) for Z ^ ro. 



Note that the function '^{r) is different from 5'(r, 0) which 
is the flux function profile on the disk surface. Hence, 
r^d<j)/dz = rBfj,{r,z)/ Bz{r,z). We take for simplicity that 
Vz = dZ/dt = const. We determine Vz in §3. In this case, 
-ff(*) = [r'^Cl{^)/V^,]Bz, where the right-hand side can be 
written as a function of \1/ and d"^ /dr. This relation allows 
the closure of equation (3) which can now be written as 

d^-^ A- lid* A (d'^V'^dn_ 

where A = [r^/ cf[{c/Vzf - 1] = {rn/cf/iV^ - 1), and 
where r= 1/(1- l///c2)i/2. 

Solution of equation (5) is facilitated by introducing di- 

mensionless variables. We measure the radial distance in 
units of the distance to the O— point of the magnetic 
field threading the disk (equation 2). We measure the flux 
function * in units of \l/o = I'oBq/2. The fields are mea- 
sured in units of Bo which is the magnetic field strength at 
the center of the disk. Thus, B,_ = (2f)-M*/df. The disk 
rotation rate is measured in units of c^/{GM) (equa- 
tion 1). In terms dimensionless variables f, ^I/, and O, 
equation (5) is the same with overbars on these three 
variables. Note also that A = f'^U'^n'^ / {V"^ - 1), where 
7^2 = (ro/r,)2 » 1. 

We consider solutions of eqn. (5) using the approxi- 
mation to eqn. (2) of *(r, 0) = which gives Cl = 
l/(a* 4- 7?.^/^ 4''^/^). These solutions are close to those ob- 
tained using the full dependence ri(^'). We then have 
Cl = l/(a, -I- TZ'^^^'i''^/^). In this approximation there 
is a unique self-similar solution to eqn. (5) for ^ ^ 
l/7^2 where n w 7^-3/2vI>-3/4 xhis solution is 'I' = 
f'^/^/[2n{r'^ - l)]2/3, with A = 2. The dependence on r is 
the same as found in the non-relativistic case by L02. The 
dependence holds for fi = [2{T'^ - l)]^/'^ /Tl < f < f 2 = 
[2(r2 - l)]i/27^i/2/33/4^ At the inner radius fi, * = l/7^2, 
which corresponds to the streamline which passes through 
the disk at a distance r = rg. For f < fi, we assume 
4' (X f^, which corresponds to B^ = const. At the outer 
radius f2, ^ = {^)max = 1/3 which corresponds to the 
streamline which goes through the disk near the O— point 
at r = ro- Note that there is an appreciable range of radii 
if 7^^/^ ^ L The non-zero field components of the Poynt- 
ing jet are E,. = -V2{T^ -lY^^Bz, B^ = -v^FB^, and 
Bz = |(f-2/3)/[27^(^2-l)]2/3, which hold for n < r < r2. 
This electromagnetic field satisfies the radial force bal- 
ance equation, dBl/dr+ {l/r'^))d[r'^{Bl- El)]/dr = 0, 
as it should. At the outer radius of the jet at r2, there 
is a boundary layer where the axial field changes from 
Bz{r2) to zero while (minus) the toroidal field increases 
by a corresponding amount so as to give radial force bal- 
ance. This gives -r2B^{r2) = {-H)raax = 2lmaxlc = 

(2/-\/3)4'(r2)/r2 which is close the relation found in the 
non-relativistic case (L02). 

3. FORCE BALANCE AT THE HEAD OF THE JET 

In the rest frame of the central object the ax;ial mo- 
mentum flux-density is Tj^, where T^fe, (j, = 0,1,2,3) 
is the usual relativistic momentum flux-density tensor. In 
the reference frame comoving with the front of the Poynt- 
ing jet, a Lorentz transformation gives T^^ = ^"^X^zz — 
2{Vz/c)nz + {Vz/cfTmi where Tzz = {E^ + Bl - 
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BD/^h — B^{4T^ - 3)/87r, where the last equality uses 
our expressions for the fields. Further, Tq^ — Ej-B^/Att = 
S2[4r(r2 - l)i/2]/8^, and Too = {E^ + + B^^)/8n = 
i?|(4r^ — l)/87r. Combining these equations gives T^^ = 

/Stt. Since T!,^ varies with radius, we take its average be- 
tween n and r2 which gives < T^^ >= B'^/[8irn'^{r^-lf]. 
assuming n <^r2- 

In the reference frame comoving with the head of the 
jet, force balance implies that the ram pressure of the ex- 
ternal ambient medium pe^t^^V^ is equal to < T^^ >. 
This assumes Vz much larger than the sound speed in 
the external medium; however, details of the shock(s) 
at the jet front are not considered here. Thus we have 
(r2 - 1)3 = B2/[87r7^2p,,tc2]. For F » 1, 

where next = Pext/fn, with in the mean particle mass. A 
necessary condition for the validity of eqn. (6) is that the 
axial speed of the counter-propagating fast magnetosonic 
wave (in the lab frame) be larger than Vz so that the jet is 
effectively 'subsonic' This condition can be expressed as 
r2 < v'^/2, where v'^ = |B'|/(47rp')^/^ is the Alfven speed 
in the comoving frame of the jet. The above-mentioned 
force- free condition is (w^)^ S> 1. The other important 
condition is that TZ^ = (ro/rg)^ ^ 1 so that the jet energy 
is extracted mainly from the disk (rather than from the 
rotating black hole). 

The derived value of F is of the order of the Lorentz 
factors of the expansion of parsec-scale extragalactic ra- 
dio jets observed with very-long-baseline-interferometry 
(see, e.g., Zensus et al. 1998). This interpretation as- 
sumes that the radiating electrons (and/or positrons) are 



accelerated to high Lorentz factors (7 ^ lO"^) at the 
jet front and move with a bulk Lorentz factor F rela- 
tive to the observer. The luminosity of the +z Poynt- 
ing jet is Ej = c J^^ rdrErB^/2 = cBplrl/i « 2.2 x 
10''''^(So/103G)2(7^/10)(M/109MQ)2erg/s, where Mis the 
mass of the black hole. 

The region of collimated field of the Poynting jet may be 
kink unstable. The instability will lead to a helical distor- 
tion of the jet with the helix having the same twist about 
the z— axis as the axisymmetric B field. A relativistic per- 
turbation analysis is required including the displacement 
current. The evolution of the jet evidently depends on 
both the speed of propagation of the lateral displacement 
and the velocity of propagation of the "head" of the jet 
Vz- Relativistic propagation of the jet's head may act to 
limit the amplitude of helical kink distortion of the jet. 

For long time-scales, the Poynting jet is of course time- 
dependent due to the angular momentum it extracts from 
the inner disk (r < Tq) which causes ro to decrease with 
time (L02). This loss of angular momentum leads to a 
"global magnetic instability" and collapse of the inner disk 
(Lovelace et aZ.1994, 1997; L02) and a corresponding out- 
burst of energy in the jets from the two sides of the disk. 
Such outbursts may explain the flares of active galactic 
nuclei blazar sources (Romanova & Lovelace 1997; Levin- 
son 1998) and the one-time outbursts of gamma ray burst 
sources (Katz 1997). 
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